Vortex Rings in two Component Bose-Einstein Condensates 



Max A. Metlitski and Ariel R. Zliitnitsky 



m 
o 
o 

(N 



(N 



> 

in 

o 
m 
o 

B 

a 
o 
o 



Department of Physics and Astronomy, 
University of British Columbia, 
Vancouver, BC V6T IZl, Canada 

We study the structure of the vortex core in two-component Bose-Einstein conden- 
sates. We demonstrate that the order parameter may not vanish and the symmetry 
may not be restored in the core of the vortex. In this case such vortices can form 
vortex rings known as vortons in particle physics literature. In contrast with well- 
studied superfluid ^He, where similar vortex rings can be stable due to Magnus force 
only if they move, the vortex rings in two-component BECs can be stable even if they 
are at rest. This beautiful effect was first discussed by Witten in the cosmic string 
context, where it was shown that the stabilization occurs due to condensation of the 
second component of the field in the vortex core. This second condensate trapped in 
the core may carry a current along the vortex ring counteracting the effect of string 
tension that causes the loop to shrink. We speculate that such vortons may have 
been already observed in the laboratory. We also speculate that the experimental 
study of topological structures in BECs can provide a unique opportunity to study 
cosmology and astrophysics by doing laboratory experiments. 

PACS numbers: 



I. INTRODUCTION 

With the recent controlled creation of vor- 
tices 0, 13 in trapped atomic Bose-Einstein 
condensates (BECs), the study of structure 
and dynamics of topological defects in weakly 
interacting superfluids has become an active 
subject of experimental and theoretical re- 
search. 

Similar topological structures in form of 
vortex lines (strings) have been observed in 
superfluid ^He and ^He, conventional super- 
conductors and high superconductors. In 
most studied and well understood cases such 
as the ^i7e system and the conventional su- 
perconductors, the order parameter can be 
represented by a single complex scalar field, 
and the structure of the vortex is quite sim- 
ple. It can be understood in terms of sponta- 
neous violation of the U{1) symmetry, which 
represents the conservation of number of par- 
ticles (in ^ife) or electric charge (in conven- 
tional superconductors). As is well-known. 



spontaneous violation of the U{1) symme- 
try leads to a non-zero magnitude of the or- 
der parameter, which represents the super- 
fluid density (in ^He) or the Cooper pair den- 
sity (in conventional superconductors). The 
structure of the vortex core in these cases 
is trivial: the order parameter must vanish 
in the center of the core, restoring the U{1) 
symmetry there, see e.g. Q. If one makes 
a ring out of such a vortex, it would quickly 
shrink and decay due to the large string ten- 
sion. The well-known quantized vortex rings, 
which have been produced and detected in 
superfluid ^ife, were stable due to their mo- 
tion that gives rise to a Magnus force equi- 
librating the string tension. We shall not 
discuss vortex ring stabilization by Magnus 
forces in the current work. 

The subject of this work is the analysis 
of a less trivial and much more interesting 
situation when the vortex rings, dubbed vor- 
tons, can be stabilized even when they do not 
move. As we discuss below, such a stability 
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cannot occur in a system described by a sin- 
gle scalar field; rather such a stability may 
occur in more complex systems, such as the 
two-component BECs, which are considered 
in this paper. 

We should note at this point that the the- 
oretical feasibility for such a phenomenon to 
occur, was first discussed by E. Witten who 
suggested a toy model leading to existence 
of superconducting cosmic strings |5|. The 
idea of construction is simple: consider a two 
component system (such as a two component 
BEC), which is described by two complex 
scalar fields ('0, </>) with an approximate U{2) 
symmetry. If the U{2) symmetry between 
fields and (p is explicitly broken down to 
f/(l) X f/(l), the condensation might be 
energetically more favorable than conden- 
sation, and the ground state will be given by 
(ip) 7^ and (0) = 0. As is well known, such 
a system allows for existence of ip vortices 
characterized by the phase of ip field varying 
by an integer multiple of 27r as one traverses 
a contour around the vortex core. What Wit- 
ten actually has demonstrated is that if the 
approximate U{2) symmetry is broken only 
weakly, the (p field may condense inside the 
core of the ip string, "spontaneously break- 
ing" the corresponding U{1) symmetry in the 
core. In the case when the (p field is elec- 
trically charged, such a string can support 
a persistent electric current along the string 
core, hence the term "superconducting cos- 
mic strings". As will be discussed below, 
the question of existence and magnitude of (p 
field condensation in the vortex core depends 
solely on how strongly the U{2) symmetry is 
broken. 

One should remark here that the above 
idea (that the core of the string could be 
in a different phase without restoration of 
symmetry) was motivated by the grand uni- 
fied theories with a typical scale of lO^^GeV. 
However, this phenomenon was experimen- 
tally observed in a very different environment 
with a typical scale of 10~^eV. Here we are 
referring to the experimental observation of 



ferromagnetic ^He — A cores in superfluid 
^He — B vortices, for review see 0,[2|- Also, 
this phenomenon has been predicted to oc- 
cur in the SO (5) model of high-temperature 
superconductivity, where the cores of the 
conventional Abrikosov strings are in the 
antiferromagnetic state 0, 0|. There has 
been recent experimental evidence that sug- 
gests this theoretical picture may be correct 
3jiJJ0 EH- It has been argued 
ISl Il6l . Il7| that this phenomenon may also 



occur in the so-called color superconducting 
phase of QCD that is believed to be realized 
when the baryon density is a few times larger 
than nuclear density jl8l |. 

The consequences of core-condensation 
are far-reaching and cannot be overseen at 
this point. In particular, as will be demon- 
strated below, one expected consequence of 
core-condensation is expansion of the vortex 
core as the magnitude of explicit U{2) sym- 
metry breaking decreases. Thus, in the limit 
of very weak U{2) symmetry breaking the 
vortex cores in 2-component BECs may be- 
come large enough to allow for direct exper- 
imental study of core properties. One other 
known consequence of core condensation is 
that the interactions between vortices can be 
drastically altered by the presence of non- 
trivial cores 3]. If this indeed happens, the 
standard picture of Abrikosov lattices in high 
Tc superconductors or vortex lattices in BECs 
may not be necessarily always correct. 

Another consequence of the phenomenon 
of core-condensation was discussed previ- 
ously in the context of cosmic strings. Within 
Witten's original toy model, it has been 
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27 that 
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condensation inside the core provides a way 
to stabilize a string loop against shrinking, 
see [ii^ for review. Indeed, a string loop is 
a topologically trivial object that normally 
shrinks due to string tension. However, if one 
allows for the (p field condensed in the core of 
a Ip vortex to carry a winding number along 
the string as well Noether's charge asso- 
ciated with its U{1) symmetry, then the con- 
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servation of both (the charge and the winding 
number) prevents the loop from shrinking. In 
fact, when both Noether's charge and wind- 
ing number are present, the field carries 
a non-zero angular momentum perpendicu- 
lar to the plane of the loop, which stabilizes 
the configuration. The resulting stable vor- 
tex loops, known as vortons, are called semi- 
topological defects since they are stabilized 
in part by topology and in part by energetics. 
One should comment here that this stability 
is purely mechanical in origin and to a lead- 
ing approximation is independent of whether 
the relevant symmetry is local or global. 

Up to this point, most discussion of vor- 
tons has been confined to relativistic physics 
(such as cosmic strings and high-density 
QCD), or to physics where the relevant low- 
energy effective Lagrangian has a relativis- 
tic form (such as the SO (5) model of high 
Tc superconductivity). In the present work 
we wish to investigate whether counterparts 
of semi-topological defects described above 
exist in non-relativistic field theory, such as 
the one describing a mixture of two ultra 
cold Bose condensed gases. In this paper we 
discuss the structure of BECs with two in- 
ternal levels. This is equivalent to a spin- 
1/2 fluid: the order parameter has U{2) ro- 
tational properties. We demonstrate that 
strings with condensation in the core do ex- 
ist if the U{2) symmetry is slightly broken, 
and we calculate explicitly under which con- 
ditions this occurs. We then go on to con- 
struct vortons in non-relativistic systems and 
demonstrate their stability under certain con- 
ditions. 

The application to the field of Bose gases 
seems particularly interesting in the wake of 
recent theoretical and experimental achieve- 
ments in this area. On the experimental side, 
the long sought goal of obtaining a weakly 
interacting Bose gas is now achieved in ultra- 
cold alkali metal gases such as Rb, Na and 
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^ ^ Moreover, experimentally, it 

is now possible 3^, S^] to have a system with 
two hyperfine states of ^'^Rb condensed in the 



same trap; this is precisely the kind of system 
we consider in the present work. Also, as we 
already mentioned, in the past few years ex- 
periments have confirmed that rotational vor- 
tices do indeed form in the ultra-cold Bose 
gases HI, 3 Q. Actually, the original ex- 
periment that produced a quantum vortex in 
^^Rb [ij for technical reasons made use of a 
second condensate in the core of the vortex. 
Thus, we speculate that the phenomenon of 
core condensation described in this paper 
might have been already observed in BECs. 
However, since in the original experiment the 
second field component was artificially in- 
duced in the core, one cannot immediately 
identify the phenomenon of natural core con- 
densation that we are discussing with the ex- 
perimentally observed result. Moreover, we 
note that vortex rings, which are the subject 
of the present work, may have also been ob- 
served However, an additional analysis 
is required before one can convincingly argue 
that the observed rings are precisely vortons, 
which are stable due to the condensation in 
the core rather than due to some other fac- 
tors such as nonzero velocity of the ring or 
due to the influence of the trap potential. 

On the theoretical side, the weakly in- 
teracting BECs provide a fascinating tool 
for testing theoretical ideas in various fields 
where experimental control is not possible 
(such as cosmology and astrophysics, see 
some comments in the Conclusion). Recently 
there has been a number of theoretical inves- 
tigations into topological defects such as do- 
main walls, vortices and skyrmions in Bose- 
Einstein condensates |36|]-||48j]. Most of these 
papers, however, are focused on the case 
when only a single field is present so that the 
effects discussed in these paper do not occur 
at all. In other papers 1.20, two particle 
species are present in the ground state. In 
this case, the effect that is the subject of this 
paper (when the second condensate is sup- 
pressed in the bulk of the media but reap- 
pears in the core of the vortex) is not as pro- 
nounced. In a configuration known as 
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a " dark- bright vector soliton" is discussed, 
consisting of a domain wall formed by one 
condensate with the second condensate con- 
fined to the wall's center. This defect is 
just a 1-dimensional analogue of the vortex 
with a core condensate considered in this pa- 
per. In other related works ji^, li^ 42, 
a nontrivial vortex ring in a two-component 
BEC is interpreted as the skyrmion, a topo- 
logically stable solution originally introduced 
in particle physics literature [i^. While the 
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analysis!^, Ha, |42l, |4^] is mainly numerical, 
and therefore, it is quite difficult to make 
a precise correspondence with our analytical 
approach, we have a strong feeling that con- 
figurations described in H Q are 
more similar to non-topological vortons con- 
sidered in the present work rather than to 
topological skyrmions as claimed, see argu- 
ments in js^l- The advantage of the analyt- 
ical approach advocated in the present work 
is the clear understanding of the structure of 
the solution as well as the nature of the sta- 
bility of the configuration. Numerical calcu- 
lations are more useful when the infiuence of 
the trapping potential (and other complica- 
tions, which are always present in real exper- 
iments) cannot be neglected and should be 
taken into account. 

This paper is organized as follows. In 
Section II we review the properties of the 
conventional U{1) strings and construct ex- 
plicit numeric and variational solutions for 
the ip string. In Section III, we determine 



the precise conditions under which string 
core condensation occurs and construct nu- 
merical solutions for and ip fields in this 
regime. In Section IV, we describe vortons 
in non-relativistic physics and demonstrate 
their classical stability. 



II. THE LAGRANGIAN AND THE 
CONVENTIONAL U{1) STRINGS 

We wish to consider a non-relativistic sys- 
tem with two types of particles such that two 
different species are very similar to each other 
and the number of particles of each type is 
conserved, i.e. the exact symmetry of the 
system is U{1) x f/(l). We also assume that 
the s wave scattering lengths between like 
and unlike species are numerically close to 
each other, such that the relevant symmetry 
is even larger, as will be explained below. Fi- 
nally, we assume here that the trapping po- 
tentials are sufficiently wide so that in what 
follows we neglect all the boundary effects. 
The last simplification is mainly motivated 
by our interest in effects, which may occur 
in very large systems (such as the systems, 
which appear in cosmological and astrophys- 
ical context) rather than in finite volume sys- 
tems where the trapping potential plays an 
essential role. Thus, throughout this paper, 
we consider the following Lagrangian describ- 
ing a system of two weakly interacting ultra- 
cold Bose-gases of same mass: 



C = zh<l>^do<l> |V$P + /ii| 

2m 



+ /^2|0P - -gn\ 



where $ = {ip,(p) is a doublet of complex 
scalar fields^. The chemical potentials of 



^ The original dimensionfuU fields and space-time 
variables appearing in Q carry a sign " tilde" . 



Throughout this paper we will be mostly using 
a dimensionless notation for the fields and space- 
time variables, which will be introduced below; 
these will carry no sign " tilde" . 
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particles {ip, (p) are (/ii,/i2) and coupling 
constants gij are proportional to scattering 
lengths aij via: gij = ^aij. By varying 

the action S = J CcPxdi with respect to ip,(f), 
one can derive the familiar Gross-Pitaevskii 
equations. 



If {il',4>) denote two hyperfine states of 
^^Rb, the scattering lengths an, ai2, differ 
by only about 3%. We neglect this difference 
in what follows, and assume a^j = a, g^j = g, 
such that the Lagrangian can be written 
in a U{2) notation as follows. 



C = zh^^do^ |V<lp- 

2m 

In the present work we assume /ii > /i2, and 
hence the ip field will condense, while the (p 
field will remain uncondensed in the ground 
state. 

It proves to be very convenient to intro- 
duce a dimensionless notation for fields, as 
well as for spacial (and time) coordinates. In 
the dimensionless notation all distances are 
measured in units of the correlation (healing) 
length ^ = {h^ /2m^iy/'^ , all frequencies are 
measured in units H\/h, while the absolute 
values of the fields \%p\'^ are measured in 
units of particle density n = fj.i/g, see Ap- 
pendix for precise correspondence. Through- 
out this paper, we use the sign "tilde" to 
distinguish dimensionless (without tilde) and 
original, dimensionfuU (with tilde) variables. 
Although we will be primarily using the di- 
mensionless notation, we will sometimes re- 
store all physical units in key formulas; we 
hope this will not confuse the reader. With 
all these remarks in mind, our starting point 
is the following dimensionless Lagrangian, 

C = t<^^do<^ -n, ^ = {tp,(P). (3) 

Here the Hamiltonian density T-C{ip,(p) has 
the form, 

n = I V<l>p + - 1)' + (4) 

where all derivatives are with respect to the 
dimensionless coordinates introduced in the 
Appendix. This Lagrangian is exactly equiv- 
alent to eq.©. It now is apparent that the 
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only relevant parameter entering the Hamil- 
tonian is the difference in chemical potentials 
of ip and (p fields, which is parameterized by 
a dimensionless factor Sm^ = {ni — /i2)//^i- 
In expression Q, $ is a doublet of complex 
scalar fields, and the Lagrangian possesses an 
exact U{1) X f/(l) symmetry with respect to 
independent phase rotations of the fields, 

f/(l) : ^ e'^^'ip, U{1) : (p e'"^(p. 

The corresponding conservation laws are 
those of the number of particles of each 
specie, Ni = J cPxip'^ip and N2 = J d^xcp^cp. 
Moreover, the Lagrangian Q possesses an 
approximate U{2) symmetry, which is explic- 
itly broken by the (5m^|0p term. In the limit 
6m the Lagrangian acquires an exact 
U{2) symmetry, 

with 0"" being the Pauli a— matrices. All 
of these symmetries are global, although lo- 
cal U{1) symmetries could also be consid- 
ered without changing the results of this 
work qualitatively. Without loss of general- 
ity, throughout this paper, we consider the 
case 5m? > 0. Then by minimizing the po- 
tential energy, we find that the vacuum ex- 
pectation value for ip field is non-zero, which 
we choose to be {ip) = 1, while {(p) = 0. Thus, 
the U{1) symmetry is spontaneously broken, 
while the U{1) symmetry remains an exact 
symmetry. Hence the vacuum manifold of 
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the system is a circle with points char- 
acterized by the phase of the field ip- As we 
know, spontaneous breaking of a U{1) sym- 
metry leads to formation of topological de- 
fects known as strings. Therefore, one should 
expect the existence of these objects for ar- 
bitrary non-zero values of parameter Srn^. 

Let us review here the elementary proper- 
ties of U{1) strings, which will be used later 
in the text. Until further notice, we will con- 
sider all the fields to be uniform in the 3rd 
(string) direction and, thus, constrain our at- 
tention to a plane crossing the string at a 
right angle. A string is a time-independent 
configuration, which exhibits a non-trivial 
mapping of the plane boundary into the 
vacuum manifold. In particular, for the U{1) 
string, as one traverses the circle out at in- 
finity in real space, the phase of the field 
makes an integer number of windings / in the 
vacuum manifold S^. Explicitly, 



2,71 J r 



(5) 



where F is some large circle in the plane cross- 
ing the string. Such configurations are ab- 
solutely stable since it is impossible to con- 
tinuously deform configurations with differ- 
ent values of topological charge / into each 
other without leaving the vacuum manifold. 
The first important feature of this configu- 
ration is that if / 7^ 0, if) must vanish at 
some point in order to be properly defined, 
see e.g.jl^. This implies that ip leaves the 
vacuum manifold somewhere and the config- 
uration certainly possesses a non-zero energy. 
Our next remark is that if the U{2) symmetry 
between and fields was exact (i.e. 5m? 
= 0) then the vacuum manifold would have 
been given by I'I'I = 1, rather than = 1. 
As is well known, there are no topologically 
stable strings possible in this C3jSG clS there 
is always a direction in the configurational 
space through which any windings of the 
field can be continuously removed. An im- 
portant consequence of the above formal fact 
is that in two component BECs the only rea- 
son for the existence of vortices is an explicit 



violation of the U{2) symmetry. Without this 
violation, the strings would immediately de- 
cay. Thus, we expect the typical character- 
istics of strings such as the core size to be 
ultimately related to the magnitude of viola- 
tion of the f/(2) symmetry, see the detailed 
discussion in the next section. 

Let us now construct explicit variational 
and numerical string configurations, which 
will be used later in the text. Formally, we 
wish to find a solution to the following equa- 
tions of motion, 

V^tP{x) = {\ij{x)\^ -l)^{x), (f) = 0. (6) 

such that the ip field satisfies the string 
boundary conditions, i.e. it carries a wind- 
ing number /. Here we have assumed that in 
the lowest energy state possessing the string 
boundary conditions, the field nowhere 
leaves its vacuum expectation value (we will 
see in the next section that this need not be 
the case). The standard cyhndrically sym- 
metric ansatz corresponding to the string 
configuration is. 



/(r 



ip{x) = 
0)^0, 



/(r)e^'^ 



oo 



(7) 



where (r, ip) are the standard polar coordi- 
nates, and /(r) is real. It is instructive to 
estimate the string tension (energy per unit 
vortex length) /i before presenting any nu- 
merical or variational calculations. By defi- 
nition, from eq. Q, the string tension is. 



/i 



d^xH 



d'x{\v^\'+lm'-ir 



(8) 
(9) 



For large r, the leading contribution to Ti 
comes from the gradient part associated with 
the phase variation of the ijj field, Ti. ~ I'^/r'^, 
and is not sensitive to the specifics of interac- 
tions. Now, performing the integration over a 
disc of radius A, we estimate fi ~ 27r/^ log(A), 
where A is a typical distance between vor- 
tices. In physical units (see Appendix) the 
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string tension is given by 

A 

f,:^ n-2nlHog{-), (10) 
2m t, 

where n = \ilj{r = oo)p is the density of the 
BEC measured far away from the vortex core. 
From this simple estimate it is clear that vor- 
tices with winding number / > 1 are not sta- 
ble with respect to decay into / vortices of 
unit winding number. Thus, throughout this 
paper we focus on simple vortices with 1 = 1. 

Now we want to construct an explicit vari- 
ational solution to a vortex with / = 1 to be 
used later in the text. Following J2l\ we take 
a variational ansatz that satisfies our bound- 
ary conditions: 

/(r) = 1 - e-^' (11) 

where is a free variational parameter. Sub- 
stituting this into the energy density (jH), in- 
tegrating over 2-dimensional space and min- 
imizing with respect to (3 we find: 

'' = W' ^=l^«-2K">s('''^) + i) (12) 

where we restored the physical dimensional 
units. Thus, we notice that in the absence 
of the second condensate (0 = 0), the radius 
of the string core is of order ^ - the typical 
correlation length in the theory. 

In order to test our variational solution 
(fTT|) we have also solved the exact eq. © with 
appropriate boundary conditions ((Tj), 

r + V-^/-(f -i)/ = o 

/(r-.0)-^0, /(r^cx))^l (13) 

numerically. We find an excellent agreement 
between our variational and numerical solu- 
tions. This gives us confidence that one can 
use either of these solutions for a more com- 
plicated problem when the configuration is 
unstable with respect to formation of a (/> field 
condensate in the core of the ip vortex, which 
is the subject of the next section. 



III. CONDENSATION IN THE 
STRING CORE 

As we know from the previous section, if 
the parameter Sm^ is zero, the Lagrangian is 
invariant under the symmetry group SU{2) x 
U{1) U{1) (broken down to U{1)). From 
topological arguments we know that such a 
Lagrangian does not allow for existence of 
vortices since the vacuum manifold is a 3- 
sphere and, therefore, does not have non- 
contractible loops. In the opposite limit, 
when is relatively large, the residual 

symmetry group is U{1) x f/(l) — > U{1) and 
the vacuum manifold is a circle 5*^, lead- 
ing to formation of stable vortex solutions. 
From these two limiting cases, it is clear 
that there should be some intermediate re- 
gion that somehow interpolates ( ClS db func- 
tion of 5m^) between the two cases. Here is 
a possible scenario, which as we shall see in 
a moment turns out to be correct. For large 
values of nothing interesting happens: 
conventional vortices of ip field do exist, the 
symmetry is restored in the string core, the 
core size is fixed by the standard correlation 
length ^ and the field vanishes everywhere 
in space. However, when Sm^ starts to de- 
crease, at some finite magnitude of (5m^, an 
instability arises through the condensation of 
field inside the vortex core. As the magni- 
tude of further decreases, the strength 
of the condensate in the core grows and 
the size of the core, where both and ip are 
far away from their vacuum expectation val- 
ues, becomes larger and larger. Finally, at 
(5m^ = the core of the string fills the entire 
space, in which case the meaning of the string 
is completely lost, and we are left with the 
situation when the U (2) symmetry of the La- 
grangian is exact: no stable strings are possi- 
ble as we know from topological arguments. 

Actually, it is of little surprise that there 
is a tendency for the field to condense in 
the vortex core. Indeed, if < < 1, the 
potential energy of the system has a global 
minimum at ^/^ 7^ 0, = 0, and a local min- 
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imum at = 0, 7^ 0. Thus, the vacuum 
structure of the theory is given by condensa- 
tion of ip rather than 0. Yet recall that in the 
core of the vortex ip must vanish, and thus the 
system is driven in the core towards the local 
minimum where the field condenses. How- 
ever, the energetic benefit of condensing in 
the core competes with the tensional energy 
needed for to relax back to its vacuum ex- 
pectation value at infinity. As was hinted in 
the previous section, the question of whether 
indeed condenses in the core is determined 
by the magnitude of [/(2) symmetry breaking 
parameter 5m?. 

After the above presentation of some intu- 
itive arguments, we proceed with the quan- 
titative discussion. Let us take the con- 
ventional string configuration described in 
the previous section by ^string = ii'ix) = 
f{r)e^'^, = 0) and formulate the question 
of whether this configuration is stable with 
respect to any small fluctuations. We must 
stress that for 5rn? > topology guaran- 



tees that small fluctuations cannot destroy 
the string, in the sense of removing the wind- 
ings of ip field. However, topology makes ab- 
solutely no predictions regarding the internal 
structure of the string - in particular, a pri- 
ori, there is no reason for = to be a stable 
solution to equations of motion in the back- 
ground of the ip string. 

We go along the standard procedure and 
expand the energy per unit length of the 
system (jH) in the ^string background to 
quadratic order in -0 and modes: 

E{ll) = iJstring + 5^^, 0) ^ /i + SE, (14) 

with /i given by (jH)). We know that the ■0 
string itself is a stable configuration so the 
6ip modes cannot possess negative eigenval- 
ues, which would correspond to an instabil- 
ity. Therefore, we concentrate only on "dan- 
gerous modes" related to fluctuations, in 
which case 6E is given by: 



6E = J d'xim' + {f{r f - l)|0p + W|0p + ^101^) 
O = -V^ + {f{rf - 1) 



Sx(p\d + 5m^)(p (15) 



where /(r) is the solution of eq. lfT^ . and we 
keep only the second order terms in assum- 
ing that the fluctuations are small. If 5E is 
a positive quantity, then the ^string is an ab- 
solutely stable configuration and modes do 
not destroy the solution discussed in the pre- 
vious section. On the other hand, if 5E is 
negative, this means that the corresponding 
is a direction of instability in the configu- 
rational space and the structure of the stable 
vortex differs from the simple ^string solution. 
Thus, the question of whether a negative 
mode exists can be answered by determin- 
ing whether the Hermitian operator O + 5m? 
possesses a negative eigenvalue. So we have 
to solve the eigenvalue problem 00 = e0, or 



more explicitly: 



l)]0(x)=e0(x). (16) 



Notice that this is just the 2-dimensional 
Schrodinger equation in a cylindrically sym- 
metric potential V{r) = f{rY — 1, which is 
everywhere negative and tends to at infin- 
ity. It is well known that such a potential 
always possesses a bound state with e < 
"s^ j. Therefore, if e -|- dm^ < then a neg- 



ative mode exists and the string obtains a 
non-trivial core. This implies that the 
f/(l) symmetry gets "spontaneously broken" 
in the core. The above formal fact has very 
important consequences such as the existence 
of a Goldstone mode capable of propagating 
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along the core of the string. Otherwise, if 
e + 5m^ > 0, the operator O + Sm"^ is pos- 
itive definite and no core condensation oc- 
curs, in which case the conventional string 
solution ^string discussed in the previous sec- 
tion is valid, and both U{1) symmetries are 
restored in the center of the core. 

We want to determine precisely the criti- 
cal value = |e| below which the vor- 
tices obtain cores. We approach this prob- 
lem numerically. Since V{x) is cylindrically 
symmetric, we assume that the lowest en- 
ergy eigenfunction 0(x) has no angular de- 
pendence. Then the operator O reduces to. 

We discretize the operator O using our nu- 
merical solution for /(r) computed in the pre- 
vious section. We then compute the lowest 
eigenvalue of the discretized operator. This 
turns out to be e ~ —0.25, which implies that 
the critical value for Sm"^ where the transi- 
tion occurs is Srn^ = (/ii — ^2)/f^i — 0.25. 
It is quite obvious that the condensation of 
the (j) field becomes even more pronounced 
when 5m^ becomes smaller. Therefore, the 
condensation of the (p field in the core of the 
ip string occurs for < ^ 0.25. In other 
words, for this range of the U{1) sym- 
metry is spontaneously broken in the core of 
the string. 

One should remark here that for Sm^ close 
to its critical value Sml ~ 0.25, the conden- 
sation in the core is small, the fourth order 
terms in ()15p can be neglected and the eigen- 
function of O, which is well localized in the 
vortex core, is a good approximation to 0(r), 
while the term ll^l"^ in eq.(fT3j) can be used 
in a variational calculation to determine the 
normalization of this eigenfunction. How- 
ever, when becomes sufficiently smaller 
than Sm^, the condensation in the core be- 
comes large, and the shape of profile function 
0(r) is infiuenced by the fourth order terms. 
Moreover, when Sm"^ becomes very close to 0, 
the magnitudes of and ip fields in the core 



become comparable (as expected from sym- 
metry considerations) and the field can no 
longer be treated weak perturbation to 
the background ip field. In this regime, '0(a;) 
and 0(x) can be found through a full varia- 
tional calculation or by solving a coupled sys- 
tem of boundary value problems for ip and 
fields. We want to obtain explicit solutions 
for the fields ip, since they will tell us how 
the vortex evolves as we approach the topo- 
logically unstable limit ^ by varying 
the external parameters /ii,/i2- It will turn 
out that by answering this apparently theo- 
retical question we will find a way to test the 
system for core condensation experimentally 
by probing the internal structure of the vor- 
tex. 

In this paper we adopt the numerical ap- 
proach to treat the above problem. The full 
Euler-Lagrange equations for fields {ip, 0) are: 

V'lP - (|0p + |^|2- 1)^ = (18) 
- (|0P + I^P - 1)0 - Sm^(l) = 

Assuming the simplest solution, ip{x) = 
/(r)e*'^, 0(x) = 0(r), equations (fTH|) reduce 
to: 

r + If - ^/ - (0' + - 1)/ = (19) 

0" + -0' - (02 + /2 - 1 + 5m2)0 = (20) 
r 

The boundary conditions on / and are 
/(O) = 0, /(oo) = 1, 0'(O) = 0, 0(oo) = 0. 
The condition 0'(O) = is necessary for 
V20 at the origin to be finite. Before we 
demonstrate the explicit numerical solutions 
to equations P9I20|) . it is instructive to calcu- 
late the asymptotic behavior of 0(r) at large 
r, which provides us with a qualitative pic- 
ture of the core size. In this case, 1 — /(r)^ 
as well as 0(r)2 are small so that eq. (|20|) can 
be linearized to obtain: 

0" + ^0' - W0 = 0. (21) 

Solution of eq.(j2H) is known to be: 

0(r) ~Ko((5mr) ~e-''"^ (22) 
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FIG. 1: Vortex profile function (solid) and 

core condensate profile function (p{r) (dashed) 
for different values of asymmetry parameter 
dm^ = {ni — /i2)//^i- Green, red and blue curves 
correspond to 5m? = 0.25,0.15 and 0.05 respec- 
tively. The quantities plotted are square roots 
of particle densities (in units of bulk density n) 
as a function of distance from the vortex core 
center r (in units of correlation length ^). 

where Kq is the modified Bessel function. 
The behavior of solution (j22p at large r sug- 
gests that the distance scale (in physical 
units), over which condenses is of order 
^/Sm rather than the correlation length ^ . 
Hence we expect that as 6m? —>■ the width 
of (f) condensate in the core grows. But (p 
and ip particles always try to minimize their 
common territory. Thus, we expect that the 
region of the string where IV'I 7^ 1 also has 
the length scale ^/Sm rather than the typi- 
cal correlation length ^. As we will demon- 
strate below, our numerical results confirm 
this intuitive reasoning. We conclude that as 
(5m^ — i> 0, the core of the string expands and 
eventually fills the whole space. So the string 
disappears in the limit Sm? = in accordance 
with topological arguments. 

Let us now explicitly test the above pre- 
dictions. We solve numerically the system 
fll9l20p for different values of Srr?. Fig. 1 dis- 
plays the numerical solutions /(r) and 0(r) 
for = 0.25,0.15 and 0.05. Essentially, 



Fig. 1 demonstrates how the densities of 
condensed particles ^n^{r)/n = \ip{r)\ and 

yjn^{r)/n = \4>{r)\ vary near the vortex core 
with distance measured in units of the corre- 
lation length ^. From Fig. 1 it is evident that 
as Sm? decreases, the magnitude and width of 
condensation in the core increases consid- 
erably, and the core of the ip vortex expands. 
Thus, the numerical calculations support our 
intuitive explanation given above. The ex- 
pansion of the vortex core due to core con- 
densation of the second component may have 
very important practical consequences. Typ- 
ically, the core size in EEC experiments can- 
not be directly observed because ^ < Ifim, 
which usually is below the resolution limit. 
We speculate at this point that one can ad- 
just the asymmetry of the system in such a 
way that the core condensation occurs. Fur- 
ther decreasing the asymmetry by adjusting 
/ii,/i2, one can make the core large enough 
(larger than the resolution limit) to allow di- 
rect detailed experimental measurements of 
vortex core properties, such as the core radius 
and the density of the second component 
filling the core. 



Summarizing the results of this section, 
we have demonstrated that in the case when 
the U{2) symmetry of Lagrangian ^ is bro- 
ken only weakly, the ip strings might obtain 
non-trivial cores. Inside of these cores, the 
(p field is condensed and the U{1) symmetry 
is spontaneously broken. We have found that 
the necessary and sufficient condition for core 
condensation is < 6m? ^ 0.25, which im- 
plies that the effect takes place even if the 
difference in chemical potentials of the two 
components is as large as 25%. The described 
effect becomes much more pronounced if the 
asymmetry parameter is small. When the 
asymmetry Sm'^ — >■ 0, the string core expands 
and eventually fills the whole space, destroy- 
ing the topological defect. 
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IV. VORTONS 

A. Stable String Loops 

The most important resuh of the previ- 
ous section is that under certain conditions 
a condensate will be trapped inside the 
%l) string. In this case, according to Gold- 
stone's theorem, a new degree of freedom cor- 
responding to the phase of the field be- 
comes a massless Goldstone boson, which can 
travel along the core of the string with arbi- 
trarily low energy cost. 

Now, we follow Witten 0, who has con- 
structed his superconducting string by mak- 
ing a loop out of the vortex with a second 
condensate in the core. Therefore, let us con- 
sider bending the if) string with a core into a 
loop of large radius R. In order to ignore cur- 
vature effects, we take R to be much greater 
than the width of the string core. Such a 
configuration is certainly unstable from the 
topological point of view since the change of 
argument of if) along a large circle surround- 
ing the loop is always (so the loop belongs 
to the trivial topological class). From an- 
other point of view, the energy of the loop 
is E = 27TfiR, where /i is the string tension 
computed in Section II. Thus, the energy E 
is minimized when R = 0. The fact that a 
simple string loop shrinks with time in the 
relativistic setting has been confirmed by nu- 
merical simulations 26] , and the same behav- 
ior should be expected in the non-relativistic 
case provided the loop does not move and the 
Magnus force does not equilibrate the string 
tension. 

Nevertheless, as was originally proposed 



21, \2 



for the case of cosmic strings j2C 
H, H H Q, a string loop can be stabi- 
lized by a persistent current carried by the 
phase of the core condensate along the loop. 
Here we discuss this original idea applied to 
our non-relativistic system. Indeed, suppose 
that the phase of the field varies by 27rN, 
as one traverses the string loop, where N is 
constrained to be an integer since is single 



valued. Then acts as a conserved num- 
ber: in order for to change, the field 
has to vanish somewhere on the string core, 
which is energetically unfavorable by the re- 
sults of the previous section. is said to 
be a semi-topological charge since it is con- 
served in part due to topological reasons and 
in part due to energetics. We can explicitly 
parameterize the field of this configuration 
as 0(2;) = 0(r)e*'^^, where r is the distance to 
the core center in the plane perpendicular to 
the string direction, so that 0(r) is the same 
as in the previous section, and z is the spa- 
tial variable running along the string core. 
Then the variation of the phase of along 
the string loop is 2nkR, so we can identify. 

In a non-relativistic system, there is al- 
ways present another conserved number - the 
Noether's charge associated with the U{1) 
symmetry, which is simply the total number 
of particles. This is given by. 



Q 



(23) 



As we demonstrated previously, it is unfa- 
vorable for particles to leave the string, so 
the conserved charge Q will be confined to 
the string core, being effectively trapped by 
the dynamics rather than by some external 
forces. 

Notice that the variation of the phase of 
along the string influences the energy of 
the configuration through the term |920p in 
the energy density. When integrated over 3- 
dimensional space, this term contributes the 
following amount into the energy of the sys- 
tem, 



AE 



'R^ 



Q. (24) 



Thus, ignoring other contributions for a mo- 
ment (detailed analysis of this problem will 
be presented in section IV C), the energy of 
the loop becomes. 



E ^ 27r/ii? + —Q. 



(25) 
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Thus, as the charges and Q are con- 
served within the configuration, the energy 
of the vortex loop has a non-trivial mini- 
mum with respect to the loop radius at Rq = 
N'^Q/iTfi. Therefore, it is expected that the 
oop will stop shrinking at the above radius 
and a classically stable configuration known 
as vorton will form. The stability is purely 
mechanical in origin and can be thought to 
be due to the angular momentum M = QN 
carried by the vorton (see section IV D). 

The formula is quite important, so 
we want to restore all dimensional factors to 
represent the final result as follows. 



E{Ro) = 37r/ii?o, Ro 



\ 2m nfl 



(26) 



where fl is determined by expression and 
Q is the physical number of trapped parti- 
cles. The most striking feature of the vorton 
is that its size could, in principle, be arbi- 
trarily large. Practically, the size is mainly 
limited by the boundary conditions specific 
to a given experiment and by our ability to 
produce a configuration with a large winding 
number N and charge Q during the forma- 
tion period. 

It is important to note that the vortons 
described here have nothing in common with 
vortex loops (observed in the superfiuid ^He, 
see 0) stabilized by Magnus force, which is 
the interaction of a moving string with its 
environment. String rings discussed in the 
present work can be stable even if they do not 
move as a whole object. The non-relativistic 
configuration bearing some similarities to the 
one described here is the twisted vortex ring 
(vorton) discussed in the superfiuid ^He — 
B system at low temperature and pressure, 
see ^ and references therein. Although it 
was originally hoped that vortons can be- 
come stable in this system, numerical esti- 
mates have shown that it is unlikely to 
happen, and such rings have never been ob- 
served in ^He — B. Actually, vortex ring 
formation and dynamics in BECs were dis- 
cussed earherlH H H H |4l|. However, 



in most cases the analysis of ring stabiliza- 
tion (as opposed to formation) was done for 
a one-component BEC, where the stability 
(or meta-stability) is due to the combination 
of ring motion and details of trapping poten- 
tial. In contrast, our main objective is the 
analysis of vortons, which can be stable in an 
infinitely large volume without moving; this 
is the kind of system, which might be of inter- 
est for cosmology and astrophysics. A close 
analogue of a vorton is the full 3D skyrmion, 
which was recently discussed in the context 



of 2-component BECsjia, |40, lill, liSi . Simi- 
larly to a vorton, the skyrmion can be repre- 
sented as a vortex ring in one component with 
the second component carrying an integer 
quanta of circulation around the ring. Nu- 
merical simulations 0, li^ of the skyrmion 
solution in BECs resemble very much the 
analytical construction of vortons discussed 
here. We suspect that the stabihty of con- 
figurations considered in ji^ [i^ is due to 
the energetics rather than topology, which 
is claimed to be the main source of stabil- 
ity, see comment sO]- We must stress that 



these simulations were performed for small 
values of topological charge only, thus, pro- 
ducing point-like objects of small size. On the 
other hand, the size of vortons is determined 
by quantum numbers Q and and could, 
in principle, be arbitrarily large as explained 
below. 



B. Relativistic Springs and Vortons 

The detailed construction of vortons in 
non-relativistic physics will be discussed in 
the next subsection. However, we believe 
that a historical aside will be instructive 
here: we now make a detor to briefiy de- 
scribe stable vortex loops, which were first 
considered in the context of cosmic strings 
0] (for ^ foil review of the subject see 
28l|). Later on, the ideas developed in the 
above works have been applied to systems 
with high baryon density and to high 
superconductors[§]. The lessons learned from 
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these papers prove to be very useful for our 
discussion of vortons in non-relativistic sys- 
tems. 

In this subsection, we consider the phe- 
nomenon of string core condensation and vor- 
ton stability in the relativistic system defined 
by the following Lagrangian, 



(27) 



where Ti is the same as before (0)). Here we 
use the relativistic notation with h = c = 
1. The Noether's charge associated with the 
f/(l) symmetry in the relativistic case is: 



Q 



d^x{do<p<f - do(p* <p). 



(2J 



It is important to note that in the relativis- 
tic setting, the (f) condensate in the core of 
the string does not necessarily carry a t/(l) 
Noether's charge. In particular, all static so- 
lutions satisfy Q = 

The first attempts to stabilize the string 
loops were based on the introduction of a 
winding number density k = N/R of the 
core condensate 0, exactly as described in 
the previous section. These relativistic ob- 
jects, dubbed springs, were characterized by 
N ^ and Q = 0. It was soon real- 
ized, however, that such objects, in gen- 
eral, are not stable: persistent current re- 
lated to the winding number cannot stabi- 
lize the string loops against shrinking. For- 
tunately, a related mechanism which makes 
string rings stable with respect to classi- 
cal decay was found by Davis and Shellard 
m H H Q. The idea is to provide 
the configuration with a nonzero charge ()28|) 
trapped in the string core along with the 
winding number N; such a configuration car- 
ries an angular momentum QN, hence its 
name - " vorton" . According to (j28|) , the pres- 
ence of charge Q automatically implies time 
dependence of the field, which can be pa- 
rameterized as 4>{x^t) = e*'-*^^~'^*''0(^)) where 
(j){r) is the profile function similar to one dis- 
cussed previously. Naively, one could think 
that time dependence of a classical solution 



brings in additional energy into the system, 
which usually further undermines stability. 
However, as one can show, the stability is 
enforced by conservation of charge Q. In a 
sense, the time-dependent configuration be- 
comes the lowest energy state in the sec- 
tor with given quantum numbers: non-zero 
charge Q and winding number N. A similar 
time-dependent ansatz for a different prob- 
lem was discussed by Coleman in js^l where 
he introduced the so-called Q-balls: macro- 
scopically large stable objects with time de- 
pendent fields. It is currently believed that 
similarly to Q-balls, relativistic vortons are 
classically stable objects. The recent, full 



scale numerical simulations|27l| support this 
statement. In the non-relativistic models 
to which our paper is devoted, a Noether's 
charge Q is always present when the field (j) 
is condensed in the core, so the vorton stabi- 
lization mechanism is very similar to the one 
suggested by Davis and Shellard. 

C. Non-relativistic Vortons 

Let us now demonstrate in detail how the 
string loop shrinks and how the vorton sta- 
bility is achieved in a non-relativistic system. 
We are interested in configurations character- 
ized by two quantum numbers. 



2ttJ c J 



(29) 

where C is the path along the core of the 
string loop. is the topological charge, and 
Q is a U{1) Noether's charge, which is just 
the total number of type particles. Fixing 
the charges Q, N, and also for now fixing 
the loop radius R, we wish to determine the 
corresponding fields. We can adopt an ansatz 
for the field motivated previously. 



0(x,t) 



^i{kz-u)t) 



(30) 



where 0(r) is the radial profile of the core 
condensate in the plane perpendicular to the 
string direction. As already explained, k acts 
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winding number density, and k = N/ R. 
As we will see shortly, u can be identified 
with the Lagrange multiplier for conservation 
of charge Q within the configuration under 
discussion. At this point, u is still a free pa- 
rameter; it will be later eliminated by solv- 
ing the equations of motion using the ansatz 
(jnn|) and substituting the solution into the 
constraint for conservation of charge Q ((221) • 
A more precise physical meaning of uj will be 
given later in the text. 

Substituting the ansatz (jHUj) into the ac- 
tion corresponding to Lagrangian 0, we ob- 
tain: 

= J dt{uj J d^x\(f){x)\^ - J d^xUm) 

(31) 

where Ti. is the full Hamiltonian density of 
the system defined by eq. (^. The Lagrange 
multiplier nature of uj is now evident - the 
energy of the system is minimized subject to 
constraint / d'^x\(j){x)\'^ = Q. In terms of the 
core condensate profile 0(r) (see eq. IHUj) . the 
above constraint requires, 

= jd'xW)\' = ^ (32) 

The physical meaning of the parameter E2 
introduced above is quite obvious: it is the 
number of particles trapped in the core per 
unit length of the string. In physical units, 

2nR J ' ' 

where Q is the physical total number of 
trapped particles and S2 is the dimension- 
less number in eq. (jS^ . 

We want to construct an effective La- 
grangian describing the "dynamics" of the 
field 0(r) on the two-dimensional xy plane 
perpendicular to the string axis, which points 
in the z direction, see (|3(J|). The effective 
Lagrangian is obtained by substituting the 
ansatz into the Lagrangian Q, 

A// = -\v^\'-lm'-i?~^'\<i>\\ 

Sm = {Sm^ + k'^ -to), (33) 




0.05 0.1 0.15 0.2 0.25 

5m^ 



FIG. 2: Number of (p particles trapped in the 
vortex core per per unit vortex length, S2 (mea- 
sured in units of n^^) function of effective 
asymmetry parameter dm = [fii — /U2 + -2:^ 

huj)/^ii. 

where the gradient is with respect to two co- 
ordinates (x, y) perpendicular to the axis of 
the string. Hence the dynamics of the system 

on the xy plane are determined by the new 

- — -2 

effective parameter 6m defined as follows: 

2 h"^ ~ ^ 

Sm = {{ni - 1x2) + - huj)/iJ,i, (34) 
2m 

where we restored in this expression all di- 
mensional factors to stress its physical signif- 
icance. This parameter plays the same role as 
the U{2) asymmetry in section III: the 
equations of motion for profile functions iplr) 
and 0(r) imposed by the Lagrangian are 
the same as eq. (jl8p . but with replaced 

by 6m : 

- (|0|2 + |^|2_ 1)^ = (35) 

- (101^ + 1^1^- 1)0 -5^V = 0(36) 

where all derivatives are with respect to co- 
ordinates {x, y) in the plane perpendicular to 
string axis. 

Following the procedure of section III, for 

—2 

each value of 5m , we numerically solve the 
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equations of motion ()35I36|) subject to string 
boundary conditions, to determine the pro- 
file functions ip{r) and 0(r). We then use the 
solution 0(r) to find E2 (see eq. ((221)) ^^s a 

function of 6m . Fig. 2 shows how the 
particle density per unit vortex length, S2 = 

S2/(^^^), depends on the effective asymme- 

— 2 

try parameter Sm . As one can see from this 
plot, the magnitude of S2 increases drasti- 

cally when 6m decreases. This is due to the 
combination of two effects. First, the magni- 
tude of the field in the core increases when 
—2 

6m decreases, as shown in Fig. 1. Second, 
the core size itself increases with decrease in 

~2 

6m . Also, as Fig. 2 demonstrates, S2 van- 
2 

ishes when 6m > 6ml ^ 0.25, since as was 
demonstrated in section III, the condensation 
of the (j) field in the core becomes unfavorable 
for large values of U{2) asymmetry parame- 
ter. 

Now we wish to understand the physical 
meaning of the parameter u, which was orig- 
inally introduced in a formal way as a La- 
grange multiplier. We define our configura- 
tion by fixing the quantum numbers Q, N 
and the loop radius R. To possess these quan- 



tum numbers the fields must satisfy the con- 
straints ((23), which fix the values of A; = 
N/R and S2 = Q/27iR. Yet, as we have 
shown above, there is a one-to-one correspon- 
dence between S2 and the effective symmetry 
breaking parameter 6m'^ (see Fig. 2). Thus, 
the value of 6m'^ is also fixed by the exter- 
nal parameters {Q,N,R). This, in turn, im- 
plies that uj is fixed via eq. (jH^ and hence 
uj is completely determined by initial assign- 
ments of {Q, N, R). Now trading external pa- 
rameters {Q,N,R) for {Q,N,k), we specu- 
late that uj{k), implicitly determined by eq. 
((Hnj). is just a complicated dispersion relation 
describing how the Goldstone mode propa- 
gates along the string loop in the presence of 
charge Q and winding number N. 

Our next task is the analysis of energetics 
of the vortex loop, with the goal of finding 
how the ring energy depends on its radius 
R. In what follows we apply the technique 
developed inli^ to the analysis of our non- 
relativistic system. Starting from the La- 
grangian Q, we switch to a Hamiltonian for- 
malism and represent the total energy of the 
configuration in the following way. 



E = 2nfiR + Jd'x{\V(P\^ + m^-l + 6m'M^ + ^\<P\'') (37) 

= 27ifiR + 27iR J d^x[(P*{-V^ + IV^P -l + 6m^ + P)0 + 

I 



where 2'KjjLR term is the conventional con- 
tribution due to the ip field, and in the 
last step we have performed an integration 
by parts. The laplacian in the last equation 
is taken with respect to the (x, y) variables 
perpendicular to the string axis. The above 
expression can be simplified via the Euler- 
Lagrange equation for the field (j){r) (j36|) . 
Multiplying by the complex conjugate 
0*(r) and then inserting the result into (p7|). 



we find: 

E = 27r/ii? + 2nR J d^x(u\(P\'^ - 
= 27r/ii? + 27ri?(cjS2 - ^24) , (38) 

where S4 = / d'^x\(j)\'^. A similar expres- 
sion for the relativistic vortons was derived 
in[28']. As we discussed earlier, u in this ex- 
pression can be computed by inverting the 
constraint ()32|33|) producing a complicated 
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implicit function of k, which is difficult to 
analyze in general. However, for relatively 
large u, such that uj > 5m^,5im? (which are 
smaller than (5m^ ~ 0.25) one can use eq. 
dnSl) to approximate u ^ k'^ = N'^/R^. The 
physical meaning of this simplification is that 
the typical energy scale of the field excita- 
tions along the string is larger than the typ- 
ical scale of excitations perpendicular to the 
string. To further simplify our discussion, we 
would like to consider the simple case when 
the core condensation is not very large and, 
therefore, the back reaction of the field onto 
the if) vortex can be neglected, so that one can 
use the unperturbed expression (fT^ for the 
string tension n Q. Also, one can check that 
in this regime S4 <C E2, so the last term in 
equation can be dropped, and the total 
energy of the configuration becomes. 



Ar2 

E = 2nfiR + —Q, 
R^ 



(39) 



where we replaced — > fc^ = N'^/R'^. 
It is interesting to note that the equation 
(j39|) exactly coincides in form with the one 
()25|1 derived with a number of simplifications 
mainly for demonstrative purposes. Thus, as 
was originally pointed out, conservation of 
charges Q and guarantees that the loop 
energy E has a non-trivial minimum with re- 
spect to radius R at. 



Rn 



vr/i 



E{Ro) = SnfxRo (40) 



The corresponding expressions with restored 
dimensional parameters were presented ear- 
lier in (j2ni)- As we argued above, we expect 
that at this radius the string loop will become 
stabilized against further shrinking, forming 
the so-called "vorton". The winding number 
density k and the particle number density E2 
corresponding to the radius Ro are: 




Due to the importance of these expressions 
as explicitly measurable quantities, we would 
like to present them with all dimensional fac- 
tors included. 



So = ^ 



3, 2m7r/i 



n 




(42) 



(41) 



One immediate observation is that if one ne- 
glects a somewhat weak InR variation of n 
with R, see f|T^ . the local parameters such 
as k and S2 depend solely on the ratio Q/N. 
Hence, the local parameters are approxi- 
mately scale-invariant under the transforma- 
tion {Q, N) — > {aQ, aN), where a is an arbi- 
trary transformation parameter. Conversely, 
the global characteristics such as the radius 
Ro and the energy E are extensive param- 
eters and transform as {R,E) — > {aR,aE). 
Thus, one should expect the existence of sta- 
ble vortons of arbitrarily large size. This 
statement, of course, should be corrected 
when finite volume effects, such as the trap- 
ping potential, are taken into consideration. 

One can present the arguments showing 
that the back-reaction of the field onto the 
ip vortex is negligible for vortons in a large 
region of {Q, N) parameter space, where, in 
particular, the parameter u is large com- 
pared to drri^.drn?. Thus, both assumptions 
made in our discussion of energetics of vor- 
tex loops with core condensates are justified 
a-posteriori. One should note that when ei- 
ther of these two assumptions is not valid, we 
still expect that stable vortons will generally 
exist. However, the analysis of vorton stabil- 
ity becomes very involved in this the 
string tension /i and the frequency uj become 
very complicated functions of {Q,N,R), and 
no analytical control is possible. We have 
not discussed this case because our main goal 
was to demonstrate the fact of vorton exis- 
tence in a large region of the parametrical 
space rather than a complete analysis of the 
allowed region of parametrical space where 
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vortons are stable. 



V. CONCLUSION 



D. Vortons and Angular Momentum 

In this section we demonstrate that the 
non-relativistic vorton carries an angular mo- 
mentum. From Lagrangian 0, the total an- 
gular momentum of the system is: 



M = j (fx(^^^ [x X -zV]$ 



(43) 



When the $ field forms a vorton, there is 
no net circulation of if) particles and the 
only contribution to the angular momentum 
comes from the field condensed in the core. 
We can re-parameterize the ansatz (IHUj) as: 



0(x) 



(44) 



where (p, 9, z) are the polar coordinates with 
respect to the center of the vortex loop. Sub- 
stituting this into gives: 



M = Nz dPx 



QNz, 



(45) 



which is a very simple relation that can be 
easily understood on the intuitive level in 
terms of classical physics. 

Therefore, the vorton carries an angular 
momentum perpendicular to the string loop 
plane. The magnitude of the angular mo- 
mentum is proportional to classically con- 
served charges Q and N . Thus, as original! 



alJT 

ii, 



suggested in the relativistic context by 
the stability of vortons can be understood in 
terms of conservation of angular momentum. 
Let us note that due to the scaling properties 
of eqs. HHl) it is energetically favorable 
for two identical vortons, each with an angu- 
lar momentum M, to merge into one vorton 
with an angular momentum 2M, thereby in- 
creasing the vorton radius (this statement is 
subject to the assumption that that the num- 
ber of (J) particles trapped in the ring core is 
conserved in this merge). So we expect that 
a typical vorton will possess, in general, the 
largest possible size constrained only by the 
boundary conditions and/or by the proper- 
ties of the trapping potential. 



The main goal of this work is to demon- 
strate that in a non-relativistic system, such 
as a two component BEC, with a U{1) x U{1) 
symmetry and an approximate U{2) sym- 
metry, stable strings with a non-trivial core 
may exist. Inside this core the U{1) sym- 
metry is spontaneously broken allowing a 
Goldstone boson to travel along the string 
axis. We have explicitly demonstrated that 
within our model the condensation inside 
the core occurs even for relatively large val- 
ues of U{2) symmetry breaking parameter, 
(5m^ ^ 0.25. The effect described becomes 
much more pronounced when the U{2) asym- 
metry becomes even smaller. We should note 
that the phenomenon when a string core pos- 
sesses a condensate of a different field/phase 
is by no means a completely new situation in 
physics. A similar phenomenon takes place 
in ^He and in high Tc superconductors where 
it has been studied experimentally quite ex- 
tensively. Nevertheless, it would be very in- 
teresting to study this very nontrivial phe- 
nomenon experimentally in two component 
BEC systems, with the specific focus on the 
dependence of the effect on the asymmetry 
parameter 5m^. In particular, we believe that 
with decrease in asymmetry, the vortex core 
will expand above the resolution constraints, 
therefore, providing an opportunity to exper- 
imentally study the core in detail. 

In the second part of the paper we ap- 
plied the idea of f/(l) symmetry breaking in- 
side the core to construct stable string loops 
known as vortons. These vortons are stabi- 
lized by a winding number carried by the core 
condensate field around the loop as well as 
a Noether's charge associated with the U{1) 
symmetry trapped in the string core. We 
have seen that the energy of such a string loop 
will have a non-trivial minimum with respect 
to the loop radius where the loop will become 
stabilized against further shrinking. We have 
also shown that non-relativistic vortons carry 
an angular momentum, which can be seen as 
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the source of their classical stability. 

It is puzzling that although, in general, 
the vorton is expected to be a fairly com- 
mon configuration in quantum field theory, 
as far as we can tell, it has never been ob- 
served in nature. Indeed, at this point, dis- 
cussions of cosmic strings [5| in cosmology and 
vortonsfl3| in high baryon density systems 
should be considered as no more than beau- 
tiful speculations. Some systems proposed 
earlier to test the existence of vortons are 
high Tc superconductors [9'] and liquid ^He. 
However, currently no direct observations can 
be presented to support the existence of vor- 
tons in high-Tc superconductors, while in the 
^He system, estimates suggest that vortons 
cannot be stabilized 0|- Thus, the BEC 
system might provide us with the first true 
experimental opportunity to prove the exis- 
tence and stability of vorton configurations 
that have been the subject of intense theo- 
retical discussions for the past 20 years. 

It is tempting to identify the vortex rings 
observed in the BEC system [35] with vortons, 
which are the subject of this work. However, 
such an identification would be premature. 
Actually, in the original experiment, a two- 
component system was initially constructed, 
but the second component was removed be- 
fore the rings were created. Hence the ex- 
perimental results were compared by 351 to 
a simulation of a one component BEC [33], 
where vortons cannot exist, as there is no sec- 
ond component capable of condensing in the 
vortex core. If, indeed, the observed rings 
are due to one field component only then 
they may be analogous to the vortex rings 



in ife|^, where the ring stability is due to 
the Magnus force. In this case, the rings must 
move in order to equalize the string tension 
that causes them to shrink. Ring stability 
might also be due in part to the artificial 



trapping potential j40|, |41|; this situation is 



not addressed in our work. 



in j35| could be also real vortons, which are 
stable due to the condensation of the sec- 
ond field in the core. In this case it would 
be the first ever observation in nature of 
such objects. There is a "smoking gun" sig- 
nature, which distinguishes between conven- 
tional vortex rings, similar to the ones ob- 
served in ^ife, and vortons, which have never 
been previously observed. This signature is 
the angular momentum (jl^ , which is always 
present for vortons and absent for conven- 
tional vortex rings^. We do not know whether 
the angular momentum of the ring (or equiv- 
alently the circulation of the second compo- 
nent in the core of the rin g) ca n be measured 
in experiments similar to [35|- Another vor- 
ton signature is its ability to stay still for a 
relatively long period of time. This distin- 
guishes vortons from the conventional vor- 
tex rings, which shrink and decay in the ab- 
sence of motion. We do not know whether 
the rings, which typically move when they 
are formed, can be stopped to observe their 
behavior. However, if either of the two exper- 
iments described above can be performed, it 
would be the " smoking gun" evidence demon- 
strating the existence of vortons stabilized 
due to core condensation. 

Finally, if vortons can, indeed, be formed 
and studied in BEC experiments, it would 
provide a unique opportunity to study 
not only condensed matter physics, atomic 
physics and optics, but also astrophysics and 
cosmology by doing laboratory experiments. 
In particular, one would finally be able to 
test mathematical ideas pioneered by Wit- 
ten, Turner, Kibble, Davis, Shellard and 
many others. Over the last few years sev- 
eral experiments have been done to test ideas 
drawn from cosmology (see the review papers 
0,0,153,153 for further details). One can fur- 
ther speculate on a possibility to study the 
formation, evolution and stability of different 



However, granted that some population of 
the second component remained in the orig- 
inal experiment, the vortex rings observed 



^ The angular momentum of a conventional ring is 
due to the contribution of the ip field into AI in 
ea. lf^ . which is identically zero. 
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topological objects (including vortons) rele- 
vant for cosmology by doing controlled ex- 
periments with BEC systems as discussed in 
the recent workshop ;56i]. 
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APPENDIX: PHYSICAL SCALES 

In this Appendix, we reduce the La- 
grangian Q to a dimensionless form 0. In 
order to do this, we introduce scaled coordi- 
nates and fields. 



t = , X = ^X, $ = y/n^, Uq 



/il 



(^) 
^2m/ii ^ 



1/2 



n 



/il 



/il 



(46) 



We remind the reader that the original (di- n.o "tilde". In the above notation, the action 
mensionfull) variables carry a "tilde", while 5* takes the form, 
the introduced dimensionless variables carry 



dtd^x{i<l>^do<l> - |V$| 



2^1 I 



r 



6m^ 



(47) 



where all derivatives are with respect to di- 
mensionless coordinates {t, x) . This is ex- 
actly (up to a multiplicative factor) the ac- 
tion corresponding to the dimensionless La- 
grangian Q used throughout this paper. 
The physical distance corresponding to one 
unit of dimensionless variable x is just the 



correlation length ^, the physical time corre- 
sponding to one unit of dimensionless vari- 
able tisl/ujQ and the physical density corre- 
sponding to one unit of |$p is just the bulk 
BEC density n. The typical values of param- 
eters in experiments on ^'^Rh are ^ ~ 0.1/im, 
n r-j lO^^cm"^, ujq ~ 5kHz 
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